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LONGITUDE REPOSITIONING OF HIGH POWER 
COMMUNICATION SATELLITES 
by Thomas A. O’Malley 
Lewis Research Center 

SUMMARY 

Communication satellites in synchronous equatorial orbit may have to be reposi- 
tioned from one stationary longitude to another. To minimize fuel consumption, the 
longitude repositioning maneuver is best effected by changing the orbit period with a 
thrust impulse collinear with the orbit velocity vector. The change in orbit period pro- 
duces a steady longitudinal drift of the satellite. When the satellite reaches the desired 
longitude, a second thrust impulse is applied to restore the synchronous condition. Im- 
pulsive thrusting, however, may not be possible if an electric propulsion system is 
used. Due to very low thrust levels, continuous thrust over many orbits may be neces- 
sary to effect a longitude change in a reasonable length of time. The velocity increment 
AV of the repositioning maneuver is a function of the change in longitude, the time of 
transit from the initial to the final longitude, and the acceleration imparted to the satel- 
lite by the propulsion system. 

For a fixed change in longitude and time of transit, there exists a minimum accel- 
eration for which longitude repositioning can be accomplished. As the acceleration 
approaches infinity, the AV approaches a value equal to one-half the value corre- 
sponding to the minimum acceleration. The nature of this variation in AV is such that 
a major portion of this possible AV reduction is obtained for an acceleration of ap- 
proximately twice the minimum. 

For a fixed change in longitude and fixed acceleration, there exists a minimum 
time of transit for which longitude repositioning can be accomplished. The AV for a 
transit time slightly larger than the minimum is significantly less than the AV for the 
minimum transit time. 

The AV required to remove any residual eccentricity at the end of the reposition- 
ing maneuver is negligible. The Earth’s triaxiality has a negligible effect on the AV 
requirement. 



INTRODUCTION 


Communication satellites in synchronous equatorial orbit may have to be reposi- 
tioned from one stationary longitude to another. For example the next generation of 
high power communication satellites may be used for communications over sections of 
Africa, Europe, and Asia. After some time, the satellite may be shifted in longitude 
by as much as 180° for Communications over North and South America. The systems 
engineer must allow for the additional fuel required for the repositioning maneuver. 

The term ’'stationwalking” is used in referring to the changing of the stationary 
position of a synchronous equatorial satellite from one longitude to another. For mini- 
mum propellant consumption the stationwalking maneuver is most effectively accom- 
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plished by changing the orbit period with a thrust impulse collinear with the orbit ve- 
locity vector. The change in orbit period produces a steady longitudinal drift of the 
satellite. When the satellite reaches the desired longitude, a second thrust impulse is 
applied to restore the synchronous condition. The second impulse is equal in magni- 
tude but opposite in direction to the first. 

Impulsive thrusting, however, may not always be permissible. Communication 
satellites of the future may use electric propulsion for attitude control and stationkeeping 
(see ref. 1) in order to keep the total propellant weight to a practical percentage of 
spacecraft weight. To avoid duplication of hardware, the same propulsion system may 
also be used for longitude repositioning. Due to very low thrust levels, continuous 
thrust over many orbits may be necessary to effect a station change. 

This report analyzes the stationwalking maneuver including the effects of nonim- 
pulsive thrusting. Design curves are included in the report for calculating the AV re- 
quirement for the stationwalking maneuver. 


SYMBOLS 


A 

A miJ AL,T 

a ‘ 

* 

Aa 

e 

e 

max 

f 


acceleration imparted to satellite by propulsion system 

minimum acceleration for a given change in longitude and time of transit 

semimajor axis of orbit 

variation of semimajor axis from semimajor axis of synchronous orbit 
eccentricity of orbit 

maximum eccentricity of orbit during nonimpulsive thrusting 
true anomaly 
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L 

AL 

L 1 

L 2 

Ar 

T 

T I 
min 1 AL, A 

t 

h 

t 2 

V 0 

AV 

AV’ 

AV e 

A ^maxl AL, A 
AV maxl AL,T 
A ^mhJ AL, A 

AV miJ AL, T 
y 

8 

AS 

*8 

AX 


satellite longitude 

change in longitude effected by repositioning maneuver 

satellite longitude before repositioning maneuver 

satellite longitude after repositioning maneuver 

variation of orbit geocentric radius from geocentric radius of syn- 
chronous orbit 

time of transit from initial longitude to final longitude 

minimum time of transit for a given change in longitude and acceleration 

time 

time duration of one thrusting period 

time duration of coast period 

velocity of satellite in circular synchronous orbit 

velocity increment required for repositioning maneuver 

velocity increment required for repositioning maneuver with triaxiality 
effects accounted for 

velocity increment required for removing residual eccentricity 

maximum velocity increment for a given change in longitude and accel- 
eration 

maximum velocity increment for a given change in longitude and time 
of transit 

minimum velocity increment for a given change in longitude and accel- 
eration 

minimum velocity increment for a given change in longitude and time of 
transit 

satellite longitude relative to nearest minor axis of Earth’s equatorial 
section 

average angular velocity of osculating orbit 

variation of average angular velocity from angular velocity of circular 
synchronous orbit 

angular velocity of circular synchronous orbit 
variation of satellite longitude from reference longitude 
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M Earth’s gravitational constant 

c p angular distance of satellite from x-axis 

to longitude of perigee measured from x-axis 

Subscript: 

0 initial condition 


STATIONWALKSNG MANEUVERS 
Impulsive Thrusting 

A synchronous satellite has an average angular velocity of 360° per sidereal day. 
To change the satellite position from longitude to longitude Lg, the angular ve- 
locity, and hence the orbital period, must be modified. For example, suppose that the 
change in longitude AL is 30° and the time of transit T from to Lg is 10 days, 
if the adjustment in orbital period is done impulsively, then the adjusted period must be 
such as to cause the satellite to drift an average of 3° per day. So the adjusted average 
angular velocity is either 363° per day or 357° per day, depending on whether the satel- 
lite drifts eastward or westward. After 10 days, when the satellite is at longitude Lg, 
the orbital period must be restored to the synchronous condition. 

The orbital period P is given by 


P = 27T 

where a is the semimajor axis and ju. is the Earth's gravitational constant. Thus the 
only way to change the orbit period is to change the semimajor axis. To produce an 
eastward drift, P must be made smaller and a must be decreased. For westward 
drift, P must be made larger and a must be increased. 

The most effective method of changing the semimajor axis is to thrust eastward or 
westward. Assume that Lg is to the east of L^. For impulsive thrusting, the first 
impulse must be directed westward in order to decrease the semimajor axis (see 
fig. 1(a)). An apogee is created at the point of thrust application. After the satellite 
has drifted to the desired longitude, a second thrust impulse is applied as shown in 
figure 1(b). The second impulse, equal in magnitude but opposite in direction to the 
first, restores the synchronous condition. The second impulse must be applied at the 
apogee of the stationwalking orbit in order to remove the eccentricity. For the case of 
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Lg to the west of Lp the direction of the thrust vectors in figure 1 would be reversed. 

Figure 2 shows orbital velocity and satellite longitude as functions of time for the 
case of Lg to the east of L^. The solid portions of the curves represent the average 
velocity and average satellite longitude. The dashed portions represent the sinusoidal 
oscillations due to the eccentricity of the stationwalking orbit. Although the first im- 
pulse reduces the orbital velocity, the average velocity of the stationwalking orbit is 
greater than the velocity of the synchronous orbit. 


Nonimpulsive Thrusting 


If low- thrust propulsion is used for stationwalking, impulsive thrusting is not 
possible. For nonimpulsive thrusting, the stationwalking maneuver is carried out in 
three phases. Let us assume that Lg is to the east of L^. Let A be the acceleration 
imparted to the satellite by the thruster. In the first phase, westward thrust is applied 
for a time t determined by AL, T, and A. The second phase is a coast phase during 
which the satellite drifts steadily eastward for a time tg. In the third phase, eastward 
thrust is applied in order to increase the semimajor axis to its synchronous value. If 
it is assumed that A is the same for both thrust periods, the third-phase thrusting 
time is again t^. Figure 3 shows orbital velocity and satellite longitude for the case of 
Lj = 0°, Lg = 170°, T = 28 days, A = 3.4x10”® g's. The oscillatory components of the 
velocity and longitude due to orbit eccentricity are not shown in figure 3. 

In order to minimize fuel consumption, the thrusting time 2t^ should be much 
smaller than the coast time tg. For very low thrust levels, however, 2t^ may be of 
the same order of magnitude as tg (see fig. 3). In the general case, for a given Lj, 
Lg, A, and T, longitude as a function of time (neglecting oscillatory components due to 
eccentricity) is given by 


L = =F(ilAIjt 2 + L t , 0<t<t., 
V 2a 


f8 A u 3 1 A | tf 

L = -)t ± + h v t, < t < t 1 + 1, 


a 


2a 


L 1 T 2 


2a 


a 


2a 


L = ±(AlAi]t 2 ±AlAl (t 2 - 2t 1 tg - 2t 2 j + Lp tj + tg < t < T 


(la) 

(lb) 

(lc) 
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When the double sign, ± or =f, is used, the upper sign is for the case of L 2 to the east 
of Lp The lower sign is for the case of Lg to the west of L-. Notice that L(t) is 
quadratic during the thrust periods and linear during the coast period. 

Equations (1) along with the equations of the following sections are derived in ap- 
pendix A. The following assumptions were made in deriving these equations: 

(1) The acceleration level A is constant through both thrusting periods. 

(2) Errors in the magnitude and direction of thrust and in the time- off and time- on 
of the thrusters are neglected. 

(3) The effect that the Earth's triaxiality has on satellite longitude is neglected. 
Appendix B presents a discussion of how the Earth’s triaxiality affects satellite 

longitude. The analysis shows that for T less than 100 days, triaxiality has a negli- 
gible effect on the repositioning maneuver requirements. 


STATIONWALKING REQUIREMENTS 
AV Requirement 

If AL, T, and A are known parameters, then one can calculate the thrusting time 
required for the stationwalking maneuver. From appendix A, the thrusting time is 
given by 


2t l- 

T -^/T 2 - — 

ALl 

V 3 

A 1 

The coast time, by the definitions of 

T and tp is 



t 2 = T - 2^ 



The velocity increment AV needed to carry out the stationwalking maneuver is 


AV = 2t x A = 



(3) 


In system design, two cases often occur. In one. case AL and T are fixed param- 
eters, and the system designer is interested in the functional relation of AV and A. 

In the second case AL and A are the fixed parameters, and the functional relation 
of AV and T is desired. These two cases are analyzed in the following sections. 
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Fixed AL and T 


If AL and T are fixed parameters, the acceleration level A must be greater 
than a minimum value A ta | AL T which corresponds to the case of continuously 
thrusting from Lj to L g ; that is, 2t^ = T. The minimum acceleration for fixed AL 
and T is given by (see appendix A) 


min' AL, T 


4a 1 AL 
3T 2 


( 4 ) 


Figure 4 is a plot of A mi J AL T as a function of AL/T 2 . By using equation (3), it can 
be shown that the partial derivative of AV with respect to A is less than zero, imply- 
ing that the AV decreases as A increases when AL and T are fixed. Equivalently, 
the AV has a maximum value AV max | AL T when A = A m j n | AT i T and a minimum 
value A V min | AL T when A = °°. From appendix A, AV min J AT T is one- half of 

AV maxl AL,T or 


AV 


mm 


- 1 AV 

AL, T ~ 2 max AL, T 


2a 1 AL 
3T 


(5) 


Figure 5 is a plot of AV max | AL T as a function of AL/T, 

The functional relation of AV and A for a fixed AL and T is given in normal- 
ized form by (see appendix A) 


AV _ A 
AV maxlAL,T A miJ AL,T 

Figure 6 is a plot of AV/AV max | AL T as a function of A /A min j AL By differen- 
tiating equation (6), it can be shown that the derivative of AV/AV max | AL T evaluated 
at A/A . I AT „ = 1 is infinite, indicating that a significant saving in AV is achieved 

HGIXT1 ^ X l 

when A is only slightly larger than A min l AL T - A major portion of the possible re- 
duction in AV is achieved when A is approximately twice A mi J AL T - Figures 4, 

5, and 6 can be used to calculate AV for given values of AL, T, and A, Figures 4 
and 5 are used to find A min | ALjT and AV max | AL>T . Knowing AV max | AL>T , 

A . | AT rp, and A, the AV can be found from figure 8. 


1 ", 


1 1 _ min 1 AL, T 


(6) 
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Fixed AL and A 


If AL and A are fixed parameters, the time of transit T must be greater than 
a minimum value T min | AL A which corresponds to the case of continuously thrusting 
from to Lg. From appendix A, the minimum transit time for fixed AL and A is 
given by 


min' AL,A 


y 4al AL 
3A 


( 7 ) 


Figure 7 is a plot of T min | AL ^ as a function of AL/A. By using equation (3), it 
can be shown that the partial derivative of AV with respect to T is less than zero, 
implying that the AV decreases as T increases when AL and A are fixed. Equi- 
valently, the AV has a maximum value AV may | AT A when T = T min | AL A and a 
minimum value AV min | AL A when T = °°. From appendix A, AV min | AL A and 

iV maxliL,A are S lvenb y’ 


A ^min' AL, A ® 


(8a) 


AV 


max' AL, A 


V 


4aA AL 


(8b) 


Figure 8 is a plot of AV max | AL A as a function of A AL . 

The functional relation of AV and T for a fixed AL and A is given in normal- 
ized form by (see appendix A) 


AV 


A ^maxl AL, A T min' AL, A 


1 - 1 / 1 - 


' min ' AL , A\ 


(9) 


Figure 9 is a plot of AV/ AV max | AL A as a function of T/T min | AL A . The functional 
relation of AV and T for a fixed AL and A is similar to the functional relation of 
AV and A for a fixed AL and T (compare eq. (9) with eq. (6) and fig. 9 with fig. 6). 


One significant difference is that AV/AV. 


M aV / AV mj4L,T-i 38 A / A mlJ 


max' AL, A 
AL, T 


0 as T/T. 


min' AL, A 


where- 


By differentiating equation (9), it can be shown that the derivative of 


AV/AV. 


max ' AL. A 


evaluated at T/T. 


min ' AL, A 


= 1 is infinite, indicating that a signifi- 


cant saving in AV is achieved when T is slightly larger than T . 


min' AL,A“ 


Figures 
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1 , 8 , and 9 can be used to calculate AV for given values of AL, T, and A. Figures 
7 and 8 are used to find T m J ALA and Knowing 

T min' AL A’ an( * T> the ^ can l° un( * from figure 9. 

If the ratio of propellant weight to spacecraft weight is assumed small, the propel- 
lant weight as a function of AV is given by 



W AV 



( 10 ) 


where W p is the propellant weight, W is the spacecraft weight, g is the acceleration 
of gravity, and I is the specific impulse. Equation (10) can be used to calculate the 

bp 

propellant weight once the AV is known. 


Orbit Eccentricity During Thrusting 

Any residual eccentricity remaining at the end of the stationwalking maneuver must 
be removed. With the assumption that the orbit is circular when thrusting begins, ec- 
centricity during the thrusting phase is given by 


e(t) = e 


max 


9 

. ^ 

sm — 

2 


where 


6 is the angular velocity of the synchronous orbit and e is 


e max 



( 11 ) 


( 12 ) 


Equation (11), which is derived in appendix A, is plotted in figure 10. If nonimpulsive 
thrusting is used, both thrusting periods should be an integral number of days so that the 
stationwalking orbit and final synchronous orbit will be circular. If the thrusting periods 
are not an integral number of days, the residual eccentricity will not be greater than 
e max‘ AV e be the velocity increment needed to remove the residual eccentricity. 

From reference 1, AV g is bounded by 


AV < AV q 
e e. 


e a 6 
_ max e 


max 


(13) 
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Substituting equation (12) into equation (13) results in 


AV e 



( 14 ) 


For A expressed in g’s and AV e in meters per second, equation (14) becomes 

AV e(-^-)- 2 - 9xlo8 | A l(g’ s) (15) 

\sec/ 

- 5 

From equation (15), AY is insignificant for very low accelerations (less than 10 g’s). 

For larger accelerations, the upper bound for AV g may become significant. It should 
be noted, however, that by a judicious choice of thrusting times, the eccentricity at 
thrust cutoff can be made much closer to zero than to e max - 


CONCLUDING REMARKS 

The purpose of this report is to analyze the methods and requirements for longitude 
repositioning of a satellite in synchronous equatorial orbit. The method of longitude repo- 
sitioning consists of an initial thrusting phase, a coast phase, and a final thrusting phase. 
The thrust vector is directed eastward in one thrusting phase and westward in the other. 

The requirements for longitude repositioning are given in terms of the velocity in- 
crement AV of the maneuver. The AV is a function of the change in longitude, the 
time of transit from the initial to the final longitude, and the acceleration imparted to 
the satellite by the propulsion system. 

For a fixed change in longitude and time of transit, there exists a minimum accel- 
eration for which longitude repositioning can be accomplished. As the acceleration ap- 
proaches infinity the AV approaches a value equal to one-half the value corresponding 
to the minimum acceleration. The nature of this variation in AV is such that a major 
portion of this possible AV reduction is obtained for an acceleration of approximately 
twice the minimum. 

For a fixed change in longitude and fixed acceleration, there exists a minimum time 
of transit for which longitude repositioning can be accomplished. As the time of transit 
approaches infinity, the AV approaches zero. The nature of this variation in AV is 
such that a major portion of the possible AV reduction is obtained for a transit time 
of approximately twice the minimum. 
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The AV required to remove any residual eccentricity at the end of the repositioning 
maneuver is negligible. The Earth’s triaxiality has a negligible effect on the AV re- 
quirement. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 5, 1971, 

164-21. 
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APPENDIX A 


DERIVATION OF EQUATIONS 

Assume that AL, T, and A are known quantities and that t^, t 2 , and AV are to 
be calculated. First we must find the change in semimajor axis due to tangential thrust- 
ing. From reference 2, the time rate of change of the semimajor axis for a nearly cir- 
cular orbit is 


^ = 2A (Al) 

dt 9 

where 9 is the average angular velocity of the osculating orbit. The acceleration A is 
assumed positive for eastward thrust and negative for westward thrust. If changes in 
9 are assumed to be small, Aa is given by 


Aa = 


2At 

9 


(A 2) 


where t is the time of thrusting. 

We now wish to find a relation between Aa and the change in average angular ve- 
locity Ad. The average angular velocity is given by 



where is the Earth’s gravitational constant. If the variations Aa and A 9 are as- 
sumed small, we obtain 



Substituting equation (A 2) into equation (A4) gives 

. • -3At t K E\ 

A 9 = (A5) 

a 
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The stationwalking maneuver is carried out in three phases. The first and third 
phases are thrusting phases, and the second phase is a coast. Let AL be the change 

a 

in longitude during the first phase, AL^ the change in longitude during the second 
phase, and AL c the change in longitude during the third phase. We now proceed to find 
AL & , AL^, AL c as functions of t^ and tg. The change in longitude AL & is given by 

AL = C 1 A 9 dt (A6) 

a JO 

Substituting equation (A 5) into equation (A6) and carrying out the integration gives 

-3At? 

AL = i (A7) 

a 2a 

During the coasting phase, A0 is constant and is equal to 

. -3At 1 

A 9 = (A8) 

a 

The change in longitude AL b is given by 

AL b = I 2 A 6 dt (A9) 

Substituting equation (A8) into equation (A9) and carrying out the integration results in 

-3ALt 9 

AL b = ^ (A 10) 


We make the convention that the sign of A is determined by the direction of thrust dur- 
ing the first phase. During the third phase, the direction of thrust is reversed. So A 6 
during the third phase, from equations (A5) and (A8), is given by 


-3At 


'A e = 


1 . 3At 


a a 


(All) 
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where t is the time elapsed in the third phase. K it is recalled that the first and third 
stages are of the same duration, then 


AL. 




Ad dt 


(A 12) 


Substituting equation (All) into equation (A12) and carrying out the integration gives 


-3At? 3 At? 

AL = 1+ i (A13) 

a 2a 


By adding equations (A7), (A10), and (A13), we obtain 


AL = AL„ + AL, + AL 
a b c 


-3Atjt 2 3Atj 
a a 


(A14) 


Now t 2 is given by 


t 2 = T - 2tj 

Substituting equation (A15) into equation (A14) yields 



Rearranging terms results in 



a AL 
3A 


= 0 


(A15) 


(A 16) 


(A 17) 


Using the quadratic formula, we solve for the thrusting time 2tj as follows: 


2t 


1 



| 4a AL 
3A 


(A18) 


From equation (A14), it is clear that AL and A are of opposite sign. So equation (A18) 
can be written as 
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(A19) 


2t i = T --*/r ! -f|f 


The AV is then given by 


AV = 2tj A 




IT - ^ / T 2 - 4a l AL 


3 A 


(A 20) 


We now wish to find AV min | ,p, the minimum AV for a fixed AL and T. By 

taking the partial derivative of equation (A20) with respect to |a| , it can be shown that 


3(AV) 

a|A| 


< o 


(A 21) 


for all permissible values of AL, T, and A. Thus equation (A20) is a monotonically 
decreasing function of | A | , implying that the minimum AV occurs when | A | is infinite. 
By taking the limit of equation (A20) as | A | approaches infinity, AV min | T is cal- 
culated to be 


AV 


min 1 AL, T 


2a AL 
3 T 


(A 22) 


Equation (A 21) implies that | AT T , the maximum AV for a fixed AL and T, 

occurs when A is equal to its minimum permissible value A min | AL The minimum 
acceleration occurs when 2t^ = T, or when the expression under the square root sign 
in equation (A19) is zero. Solving for A min | AL T § ives 


min 1 AL, T 


4a 1 AL 
3T 2 


(A 23) 


By evaluating equation (A20) when A = A min | rp, we obtain an expression for 
AV maxl AL,T 


AV 


max 1 AL, T 


4a | AL 

3 T 


Equations (A22) and (A24) imply that 


(A 24) 
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(A 25) 


AV. 


max 1 AL, T 


= 2AV , 


min 1 AL , T 


Combining equations (A 20 ), (A23), and (A24) gives 


AV 


AV max' AL,T A min' AL,T 


1 -4 /l - AL, T 

A 


(A 26) 


We now wish to find AV,„„| AT A , the maximum AV for a fixed AL and A. By 

max 1 AL, A’ J 

taking the partial derivative of equation (A 20) with respect to T, it can be shown that 


a(AV) 

ST 


< 0 


(A 27) 


for all permissible values of AL, T, and A. Thus equation (A20) is a monotonically 
decreasing function of T, implying that the maximum AV occurs when T is equal to 
its minimum permissible value, T min | The minimum time occurs when the ex- 

pression under the square root sign in equation (A19) is zero. Solving for T m ^ n | ^ ^ 
results in 


T 


min 1 AL , A 



(A 28) 


Evaluating equation (A20) when T = T m - n | ^ we obtain the following expression for 
A ^maxl AL, A' 


AV. 


max 1 AL 




4a A AL 


(A 29) 


The minimum AV occurs when T is infinite. Taking the limit of equation (A 20) as T 
approaches infinity gives 


A ^muJ AL, A ® 


Combining equations (A20), (A28), and (A29) yields 


AV 


T 


AV 


TJX 

max AL,A min AL,A 


1 -^t /1 _i min' AL,A 


(A 30) 


(A31) 
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Equations will now be derived for finding orbit eccentricity during the time of 
thrusting. The coordinate system adopted is shown in figure 11, where the x-y system 
is an inertial reference with origin at the Earth’s center. The x-y plane is the equato- 
rial plane. The angle co is the longitude of perigee measured from the x-axis, f is the 
true anomaly, and cp is the angular distance of the satellite from the x-axis. From 
figure 11 the angles co, f, and <p are related by the equation 

f = <p - co (A32) 

The point B on the Earth’s equator is assumed to lie on the positive x-axis at t = 0. 
Thus, at time t, the angle between the line OB and the positive x-axis is 8 & t, where 
0 is the angular velocity of the Earth’s rotation (360° per sidereal day). Let AX(t) be 
the variation at time t of the satellite longitude from the longitude of point B. Then 
AA(t) is given by 


AX(t) = <p( t) - 0 e t (A 33) 

We now proceed to find eccentricity as a function of time with thrust beginning at 
t = 0. We first derive equations for d(AX)/dt, de/dt, and dco/dt. For an orbit with 
small eccentricity, d</>/dt is given to first order in e by (ref. 2) 

^ = 6 + 2e9 cos f (A 34) 

dt 

where 6 is the average angular velocity of the osculating orbit. The variation in 6 is 
given by 


A 9 = 6 - L (A35) 

V 

6 © 

Using equations (A32) to (A35) and recalling that 9 « 6 result in 

V 

= A9 + 2e0 o cos (AX + it - co) (A36) 

dt e e 

From equation (A 5), A 8 is given by 


A 9 = 


-3At 


a 
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where A is the acceleration and a is the semimajor axis. The acceleration A is 
assumed positive for eastward thrust and negative for westward thrust. Equation (A36) 
can now be written as 


d(AX) _ i3At + 2 q q cos( AX + fi t - c o) 
dt a e e 

From reference 2, de/dt and dco/dt are given to first order in e by 

de _ 2A cos 1 _ 2 ± C0S(AX + ~ m) 

dt a6 e a6 e 

do. . 2A sin f _ 2A ein ^ + V - "> 

<* • 

dt ea0. ea0 o 

6 © 


(A 37) 


(A 38) 


(A39) 


Equations (A37), (A38), (A39) from a coupled system of three first-order differen- 
tial equations. If we assume that the change in AX is small over a one-orbit period, 
then equations (A38) and (A39) form a coupled system of two equations to be solved over 
a one- orbit period 


de 

dt 


2A cos(AX + 0 g t - co) 


a0. 


do; 

dt 


2A sin(AX + fi t - w) 

ea 0 e 


(A40) 


(A41) 


In equations (A40) and (A41), AX is assumed to remain approximately equal to its initial 
value AXq, To this point, we have not specified initial conditions. We now impose the 
initial condition that eg = 0. To determine cog, consider the case where the thrust is 
in the eastward direction. At t = 0, the satellite’s angular distance from the x-axis is 
found from equation (A33) to be AXq = <p(0). An eastward tangential thrust applied to a 
circular orbit at the angular position AXq will create an instantaneous perigee at that 
position. Thus, u)q = AXq. For thrust in the westward direction, cog = AXq + jr. The 
solution to equations (A 40) and (A 41) satisfying the initial conditions is found by assuming 
co is a linear function of time 


18 



(A42) 


e(t) = 


4 1 A 



. V 

sin 

2 


6t 

w(t) = — + w Q (A43) 

2 U 


Substituting equation (A42) into equation (A37) and integrating, we obtain the following 
solution for AA(t): 


AX(t) = AA q 


+ -3A (i - cos i t ) 
2a a el 

e 


(A44) 


As an application of equations (A42) to (A44), consider first the case where t = 0 
corresponds to the beginning of thrusting. Equation (A42) implies that eccentricity is a 
half sine wave over a 24-hour period with a maximum value of 


e 

max 



(A45) 


At the end of the 24-hour period, the eccentricity is again 0. We now re-initialize the 
problem, the only change being that AXq is no longer <p( 0). From equation (A44), the 
new AXq is given by 


AA n = <p{ 0) - — (24 hr) 2 
u 2a 


<p( 0) - 


6tt 2 A 



(A 46) 


Again the solution for e is a half sine wave over the second 24- hour period. By con- 
tinuing the process, the eccentricity will be a succession of half sine waves with ampli- 
tude given by equation (A45). 


19 



APPENDIX B 


EFFECT OF TRIAXIALITY ON STATIONWALKING MANEUVERS 

An analysis of the effect of the Earth’s triaxiality on a synchronous satellite is 
given in reference 3. A brief synopsis of that analysis is presented here. The Earth’s 
equatorial cross section is approximately an ellipse whose minor axis passes through 
74.6° east longitude and 105.4° west longitude (ref. 4). These two longitudes are 
stable points. Neglecting other perturbations, a synchronous satellite placed at either 
of these longitudes will tend to stay there. If a satellite is positioned at any other lon- 
gitude, it will undergo a longitudinal oscillation about the nearest minor axis. A simul- 
taneous oscillation in orbit radius also occurs. The period of the oscillations in longi- 
tude and radius is greater than 2.2 years. 

In this appendix, the symbol y will denote longitude relative to the nearest minor 
axis of the Earth’s equatorial section. Figure 12 presents the oscillations in radius Ar 
and longitude y for two cases. The initial longitude is 45° in the first case and 
25° in the second case. The period of oscillation is 2. 7 years for = 45° and 2. 3 
years for y q = 25°. 

In figure 13, we present another Ar - y curve, except now we neglect triaxiality 
effects and consider what happens when stationwalking from longitude y^ to longitude 
yg. We assume that the thrusting time is much smaller than the coasting time. If 
T = 30 days is the given time of transit from y^ to yg, then the first thrusting period 
would produce a change Ar in the orbit radius. This maneuver corresponds to moving 
from point A to point B in figure 13. During the coast phase, we move from B to C. In 
the second thrusting period, the synchronous radius is restored and we move from C 
to D. Now assume that the time of transit is T’ = 15 days, where T' < T. Then the 
change in radius Ar* will have to be greater than Ar. The path followed is AB’C’D 
instead of ABCD. 

Figure 14 presents the Ar - y curve when the triaxiality effect is taken into ac- 
count. Again assume that the time of transit is T or T’. The triaxiality produces a 
change in radius during the coast phase. So the first change in radius, Ar^ or Ar^, 
will not be the same as the second change, Arg or Ar^. The dashed portion of the 
curve represents the path that would occur if there were no second thrusting period. 
Here it is assumed that Ar^, Ar^, Arg, Ar^ are less than 30 kilometers. For changes 
in radius greater than 30 kilometers, the triaxiality effect is not great enough to 
"capture” the satellite and it would continue drifting westward in the absence of a sec- 
ond thrusting period. 
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Assume that y^, y g, and T are given parameters. Let Ar be the necessary 
change in radius , neglecting the triaxiality effect. Let Ar^ and Arg be the first and 
second changes in radius when triaxiality is taken into account. Let AV be the velocity 
increment for the stationwalking maneuver when triaxiality is neglected. Let AV’ be 
the velocity increment when triaxiality is taken into account. Table I presents an array 


TABLE I. - DATA FOR CASE OF '= 20° AND r 2 = 10° 


Time, 

T, 

days 

Velocity 

increment 

when 

triaxiality is 
neglected, 
AV, 
m/sec 

Velocity 

increment 

when 

triaxiality is 
accounted for, 
AV T , 
m/sec 

Change 
in radius 
when 

triaxiality is 
neglected, 
Ar, 
m 

First change 
in radius 
when 

triaxiality is 
accounted for, 

Ar r 

m 

Second change 
in radius 
when 

triaxiality is 
accounted for, 

Ar 2> 

m 

26.0 

2.19 

2.18 

30 000 

29 000 

30 700 

46.2 

1.23 

1.22 

16 900 

15 200 

18 200 

66.4 

.86 

.84 

11 700 

9 300 

13 700 

86.5 

.66 

.64 

9 000 

5 700 

11 600 

106.7 

. 53 

.50 

7 300 

3 300 

10 600 

126.9 

.45 

.41 

6 100 

1 100 

10 200 


of values of T, AV, AV', Ar, Ar^, and Arg for the case y ^ = 20° and yg = 10°. The 
thrusting time is assumed to be much smaller than T. As T becomes larger, the dif- 
ference between Ar^ and Arg becomes larger. Notice, however, that AV is only 
slightly larger than AV’ in all cases, the reason being that the sum of Ar^ and Arg 
is nearly equal to 2 Ar in all cases. 

One would suspect that AV - AV* is greatest when y ^ = -y g. In this case the 
stable longitude is midway between y^ and yg, and AV ! = 0 if T is equal to one-half 
the period of oscillation. But T would have to be greater than 1 year in such a case. 
For T less than 100 days, AV and AV* are nearly the same, even for the case 
= ~ y 2‘ 15 presents curves of AV and AV ! when y^ = -yg. 
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Figure 1. - Stationwalking maneuver with impulsive thrusting. 




Figure 2 . - Orbit velocity and satellite longitude as functions of time for 
impulsive thrusting maneuver. L 2 is east of Lj. 
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Figure 3. - Orbit velocity and satellite longitude as functions of time 
for nonimpulsive thrusting maneuver. L 2 is east of Lj. Oscilla- 
tory components due to eccentricity are not shown. 



Figure 4. - Minimum acceleration for a fixed change in longi- 
tude and time of transit as function of ratio of change in 
longitude to square of time of transit. 
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Ratio of velocity increment to maximum velocity increment 
for a fixed change in longitude and time of transit, 

AV/AV max |ALj Maximum velocity increment for fixed change in longitude 

r- and time of transit, AV max | AL j, m /sec 





Maximum velocity increment for fixed change in longitude Minimum time of transit for fixed change in longitude and 

and acceleration, AV max | AL#A , m/sec acceleration, T min | AL A , days 



Figure 7. - Minimum time of transit for fixed change in longitude and acceleration as a function of the 
ratio of change in longitude to acceleration. 



Figure 8. - Maximum velocity increment for fixed change in longitude and acceleration as a function of 
the product of change in longitude and acceleration. 
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Ratio of velocity increment to maximum velocity increment 
for a fixed change in longitude and acceleration, 
AV 'AV max lAL.A 


1 . 0 , 



Ql 1 1 I 1 1 1 1 t I ! I ! I I I I I I 

1 2 3 4 5 6 7 8 9 10 

Ratio of transit time to minimum transit time for a fixed change in longitude and acceleration, a 

Figure 9. - Normalized velocity increment as function of normalized transit time. Normalization is with 
respect to fixed change in longitude and acceleration. 



Figure 10. - Normalized eccentricity as function of time during 
thrusting period. e/e m;tY = |sin(iu/2)|: e maY = 4IAl/aejL 
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Figure 13. - Variation- in radius as function of longitude for two stationwa Iking maneuvers. Triaxiaiity effects are 
neglected. 


Variation in radius, Ar 



Figure 14. - Variation in radius as function of longitude for two stationwalking maneuvers. 
Triaxiaiity effects are accounted for. 





